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ABSTRACT 

High resolution MERLIN observations of a newly-discovered four-image gravitational 
lens system, B0128+437, are presented. The system was found after a careful re- 
analysis of the entire CLASS dataset. The MERLIN observations resolve four com- 
ponents in a characteristic quadruple-image configuration; the maximum image sep- 
aration is 542 mas and the total flux density is 48 mJy at 5 GHz. A best-fit lens 
model with a singular isothermal ellipsoid results in large errors in the image posi- 
tions. A significantly improved fit is obtained after the addition of a shear component, 
suggesting that the lensing system is more complex and may consist of multiple de- 
flectors. The integrated radio spectrum of the background source indicates that it 
is a GigaHertz-Peaked Spectrum (GPS) source. It may therefore be possible to re- 
solve structure within the radio images with deep VLBI observations and thus better 
constrain the lensing mass distribution. 
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1 INTRODUCTION 

Building upon the successful methodology of the Jodrell 
Bank-VLA Astrometric Survey (JVAS; Patnaik et al. 1992; 
Browne et al. 1998; Wilkinson et al. 1998) which is known 
to contain 6 examples of gravitational lensing (King et al. 
1999) within a parent sample of 2384 sources, the Cosmic 
Lens All-Sky Survey (CLASS; for example Browne 1999) 
aims to find new examples of strong gravitational lensing 
within a flux-limited sample of ~ 12 000 flat-spectrum ra- 
dio sources observed with the VLA in A configuration at 
8.4 GHz. Upon completion of CLASS the lensing frequency 
within the combined samples of JVAS and CLASS will allow 
constraints to be placed on the cosmological constant J7a. 
This has already been achieved for the JVAS sample (Helbig 
et al. 1999). Also, individual lens systems can be used to de- 
termine the Hubble parameter (Refsdal 1964) by measuring 
the time delay between individual images in a system. 

So far 18 lens systems have been found in JVAS and 



CLASSQ Prior to the final group of CLASS observations 
which are currently being followed up with MERLIN and 
the VLBA, this paper reports on the discovery of a new lens 
system, B0128+437, originally observed in the first phase 
of CLASS. The system was found after a re-analysis of the 
entire CLASS dataset. 



2 VLA RE-ANALYSIS AND MERLIN 
OBSERVATIONS 

In order to ensure uniformity over all epochs of observations, 
the entire CLASS dataset was re-flagged and re-calibrated 
together in a standard way within AIPS. The data were 
then re-mapped using an automatic script within DIFMAP 
(Shepherd 1997; Pearson et al. 1994) and Gaussian model 



* JVAS and CLASS lens survey papers can be found at URL 
tittp : //gladia. astro . rug . nl : 8000/ceres/papers/papers . html 
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Table 1. Image parameters for B0128+437 derived from the 
Gaussian components fitted to the VLA 8.4 GHz data. Im- 
age positions are offset from component A at (J2000.0) R.A. 
01 h 31 m 13!471 dec. +43°58'12'.'938. 
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'arcsec) 




Component 


East 


North 


Ssa (mjy) 


A+B 


+0.000 


+0.000 


14.8 
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+0.497 


-0.188 


3.9 


D 


+0.076 


-0.266 


5.S 



components were fitted to the components in each map. All 
sources with multiple compact model components with sep- 
arations in the range C/.'3 to 6"and flux density ratios ^10:1 
were selected. The beam size of the VLA at 8.4 GHz in A 
configuration is ~0'.'2, and a compact component is defined 
to have a Gaussian diameter (FWHM) ^170 mas. Also, only 
those systems where the total flux density in all components 
added together is Js20 mjy were selected. The vast major- 
ity of lens candidates found this way have been followed 
up during previous phases of CLASS. The small number of 
new candidates arising from the re-analysis and re-selection 
process were followed up with MERLIN at 5 GHz. 

The top-left panel of Figure | shows the VLA 8.4 GHz 
map of B0128+437 after the re-calibration of the CLASS 
data, obtained on 1994 March 05 as part of CLASS "phase 
1". The maximum component separation is 534 mas which 
is less than three times the beam size and as a result only 
three components are resolved. It can be seen a posteriori 
that this map is consistent with a four image lens system 
where the brightest component in the VLA map consists of 
an unresolved blend of two images. The total flux density of 
the system is 24.6 mjy and the flux density ratio between the 
brightest and faintest components is ~3.75:1. The details of 
the three components are given in Table pi. 

B0128+437 was observed at 5 GHz on 2000 March 09 
by MERLIN as part of 'snapshot' observations of the new 
lens candidates from the CLASS VLA data re-analysis. The 
data were edited and calibrated using the standard MER- 
LIN programs and AIPS, and then mapped to an rms noise 
level of 277 /xJy beam -1 with a beam size of ~50 mas us- 
ing DIFMAP. Four components were detected (Figure hi) 
with a maximum component separation of 542 mas. The 
brightest component in the VLA map is well resolved into 
two components in the MERLIN map with a separation of 
136 mas. The positions of all four components are consistent 
with those expected from a four-image lens system. The in- 
dividual components have been modelled with four Gaussian 
components giving a total flux density of 47.7 mjy. All four 
components are unresolved i.e. < 30 mas to within the errors 
and modelling the source with four 8 functions or four Gaus- 
sian components resulted in similar values of \ 2 - Details of 
the model components are given in Table B. Components are 
labelled (anti-clockwise from the brightest) A to D. 

The flux density ratios between components appear to 
have slightly different values at 8.4 GHz and 5 GHz. How- 
ever, components A, B, and D are contained within an area 
only two beam-widths across in the VLA 8.4 GHz map and 
significant errors in the flux density of each component can 
therefore arise when the map is deconvolved with the beam. 
Furthermore, the 3<r level at each of the two frequencies 
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Figure 1. MERLIN 5 GHz map of B0128+437 taken on 
2000 March 09 centered on (J2000.0) R.A. 01 h 31 m 13?420 dec. 
+43°58'13'.'02. The contours are at -3, 3, 6, 12, 24, and 48 times 
the rms noise of 277 fi3y beam -1 . The beam size is 58.8 mas 
X 47.7 mas in position angle 50?9. The top-left panel shows the 
VLA 8.4 GHz CLASS map taken on 1994 March 05 for compar- 
ison. Contours are at -3, 3, 6, 12, and 24 times the rms noise of 
271 /xjy beam -1 . The beam size is 0'.'246 X 0!'225 in position 
angle 14? 7. Note that the scale of the top-left panel is ~ 0.5 that 
of the main map. 



Table 2. Image parameters for B0128+437 derived from the 
Gaussian components fitted to the MERLIN 5 GHz data. Im- 
age positions are offset from component A at (J2000.0) R.A. 
01 h 31 m 13!405 dec. +43° 58' 13'.' 14. 
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-0.250 
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is ~1 mjy beam - . The combination of these two factors 
means that the differences between the flux density ratios 
at each frequency are not significant. 



3 LENS GALAXY MODELLING 

An attempt has been made to model the B0128+437 lens 
system with a singular isothermal ellipsoid mass distribu- 
tion (for example Kormann, Schneider & Bartelmann 1994). 
The four MERLIN image flux densities and positions have 
been used as input parameters, providing eleven constraints 
(eight from image positions and three from image flux ra- 
tios). The model parameters which were allowed to vary were 
the source and lens positions, the lens velocity dispersion, 
ellipticity and position angle. Thus there were seven free pa- 
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Figure 2. An attempt at modelling B0128+437 with external 
shear. The model images are shown as solid contours, whilst the 
time-delay surface is shown as dashed contours. The solid line 
shows the critical curve while the top-right sub-image shows the 
caustics where multiple images are formed. The source position 
is shown in the sub-image as a star to the left of the centre of 
the lensing galaxy, just inside the inner caustic of the system. 
The central cross-haired ellipse indicates the galaxy position and 
orientation. 



rameters, and hence four degrees of freedom. A 3 mas error 
in the relative positions of each component and a 20 per cent 
error in the flux densities to take account of source variabil- 
ity and systematic errors were allowed. A lens redshift of 
zi — 0.5, source redshift of z 3 = 1.5, and a flat universe with 
Qm = 1-0 were assumed. The resulting model suggests that 
the lens mass has an axial ratio b/a ~ 0.76. However, the dif- 
ferences in the positions of the four model images from those 
in the actual map are all in the range 30 to 100 mas which 
clearly shows that a singular isothermal ellipsoid model does 
not reproduce the observations well. This suggests that the 
system may consist of a complex deflector and may involve 
other nearby galaxies. 

Modelling with the addition of an external shear com- 
ponent resulted in image positions and flux ratios being re- 
produced well (see Figure fe|) with image positions within 
1 mas and flux ratios within a few per cent of the observed 
values; the shear has a value of ~0.23. This result also sug- 
gests that there is an additional lensing component external 
to the main lensing galaxy to the north-east or south-west. 
Optical and/or infrared imaging must now be carried out to 
determine the positions and numbers of possible deflectors 
before more accurate modelling can be achieved. No optical 
counterpart has been found at the position of the radio im- 
ages on POSS-II, placing a limit of tur > 20 on the main 
lensing galaxy. 

The image configuration in B0128+437 is similar to that 
of another CLASS lens system, B1555+375 (Marlow et al. 
1999). This system has been mapped with greater sensitiv- 
ity than B0128+437 and has a well-constrained mass model 



Frequency (Hz) 

Figure 3. The radio spectrum of B0128+437. The data have 
been obtained from the following surveys (with epoch) from 
325 MHz to 8.4 GHz: WENSS (05/93); NVSS (12/93); GB6 
(~/87); CLASS (03/00); CLASS (03/94). Error bars are shown 
at lcr. 



for the lensing galaxy. From the model for B1555+375, the 
magnification of individual images ranges from 0.5 to 8.2, 
and the combined magnification for all four images is ~20. 
This comparison is used in the next section. 



4 THE NATURE OF THE RADIO SOURCE IN 
B0128+437 

The integrated spectrum of the background source in 
B0128+437 can be obtained from the current MERLIN and 
VLA observations at 5 GHz and 8.4 GHz, and the GB6 (Gre- 
gory et al. 1996), NVSS (Condon et al. 1998), and WENSS 
(Regelink et al. 1997) surveys at 4.85 GHz, 1.4 GHz and 
325 MHz respectively. A plot of the radio spectrum from 
these data is shown in Figure H. Despite the sparsity of 
points, the spectrum clearly peaks at ~1 GHz which shows 
that the source is a GigaHertz Peaked Spectrum (GPS) 
source. 

Snellen et al. (2000) found that both the turnover fre- 
quency i^pcak and the turnover flux density .Speak of GPS 
sources correlates with their overall maximum angular size 
(component separation, not the maximum size of a single 
GPS source component) over a flux density range for the 
turnover ranging from ~40 mjy to ~5 Jy. It is estimated 
that i/peak for B0128+437 is ~1 GHz with a flux density 
•Speak ~100 mjy. However, the source is being magnified and 
lensed into four images. Assuming that the total magnifica- 
tion in B0128+437 is similar to that in B1555+375 (see Sec- 
tion Kl), the 'real' turnover flux density of B0128+437 is es- 
timated to be ~5 mjy. The observed near-linear correlation 
between the overall maximum angular size and S. peak^peak 
suggests that the intrinsic angular size of the radio source 
in B0128+437 is between 1 and 10 mas. Since the source is 
being magnified by up to ~10 times in the strongest image 
the angular size of the source could increase to between 3 
and 30 mas in this instance. The fact that the strongest im- 
age is unresolved in the MERLIN 5 GHz map also places an 
upper limit of ~30 mas on the apparent angular size of the 
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source in this image. VLBI observations at 1.4 GHz with a 
resolution of a few milliarcseconds might be able to resolve 
the images, which would place additional constraints on the 
lens model. 

The prospect of using this system to measure the Hub- 
ble parameter is unfortunately not promising. If the radio 
source is indeed a GPS source, it will be a member of one of 
the least variable classes of extragalactic radio source ( Aller, 
Aller & Hughes 1992). The GB6 survey data show that the 
4.85 GHz flux density was 47 ±5 mjy between 1986 Novem- 
ber and 1987 October. The MERLIN observations at 5 GHz 
on 2000 March give a total flux density of 47.7 ± 0.3 mjy. 
While this does not rule out the possibility of detecting vari- 
ability within the source, it is consistent with the source 
being non-variable. 



5 SMALL QUAD LENSES 

Since, in the VLA 8.4 GHz CLASS map of B0128+437 im- 
ages A and D were not resolved, as was also the case with 
B1555+375, other compact four-image lens systems may 
have been missed in the search of JVAS and CLASS if two 
of the images had a separation smaller than the beam size. 
In such cases a single extended Gaussian component would 
have been fitted to these images during the automatic map- 
ping process. The lens candidate selection criterion stipu- 
lating multiple compact components would then result in 
rejection of such a system as a lens candidate unless the 
remaining two images were detected and fitted with com- 
pact Gaussian components. In order to check that no such 
systems have been missed, a new set of selection criteria 
based on the expected image configurations of small-quad 
lens systems was devised. Five new candidates were found, 
all of which have been followed-up with MERLIN at 5 GHz. 
In each case, the 'resolved' component in the VLA 8.4 GHz 
CLASS map was shown by MERLIN to be a single extended 
component. All five candidates were rejected as lenses based 
on surface brightness and morphological arguments. It is un- 
likely that there are additional lens systems like B0128+437 
or B1555+375 within the JVAS and CLASS samples. 



6 SUMMARY AND FUTURE WORK 

A new gravitational lens system, B0128+437, has been dis- 
covered in the course of a re-analysis of the entire CLASS 
dataset. The system consists of four images of a background 
radio source which appears to be a GPS source. It might 
be possible to detect extended structure if the system is ob- 
served with VLBI. Initial attempts to model the lens mass 
with a singular isothermal ellipsoid distribution have proved 
unsuccessful, indicating that there may be multiple deflec- 
tors contributing to the lens. 

The discovery of B0128+437 has increased the number 
of known compact lens systems. Whilst there is no optical 
information to constrain the lensing mass, Turner, Ostriker 
& Gott (1984) suggest that lower separation lenses are more 
likely to be due to spiral galaxies. At present the mass distri- 
butions of these systems are not well understood at cosmo- 
logical lensing redshifts (for example Keeton & Kochanek, 
1998) making B0128+437 of particular interest. 



In order to constrain the lens mass model, further radio 
observations with MERLIN and VLBI must be made to bet- 
ter determine the image flux densities and positions. The de- 
tection of extended structure within the source will provide 
further parameters with which to constrain the lens mass 
model. Optical and/or infrared imaging to locate the lens- 
ing mass(es), and spectroscopy to find the lens and source 
redshifts are also needed for a satisfactory lens model. 
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